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5.4 Viscous effects on the instability of parallel flow 

The instability of parallel flows (in pipes, channels, boundary layers, jet wakes, plumes) is 
important to understand the transition to turbulence and has been studied by numerous 
theoreticians and experimentalists. The inviscid theory does not tell us why some flows 
without the point of inflection are unstable, as observed in pipes by Reynolds. 

It turns out tha viscosity can be destabilizing. This topic is very fully described in 
books by Lin, Drazin and Reid, and Drazin. We only give a brief sketch for the steady 
two-dimensional boundary layer on an semi infinite flat plate ( Blasius' problem) here. 

5.4.1 Role of viscosity 

By considering a control volume one wavelength long, it is possible to show that the balance 
of mechanical energy requires that 

— ^pM-^iN (5.4.1) 

where 

E = ^jJ {u'^ + v'^)dxdy (5.4.2) 

is the kinetic energyof the disturbance in the volume, 

pM = -pA J J u'v'—dxdy (5.4.3) 

is the rate of working by Reynolds stress against the shear in the basic flow, and 

nN^nJj {^)dxdy (5.4.4) 

is the rate of dissippation, with 

^1 dv' du' 

dx dy 

being the vorticity. Derivation of this formula is left as a homework exercise. 

For large /x, ^ < 0, hence viscosity is stabilizing. However for small the first term 
which is also affected by viscosity can overwhelm the second term near the wall, hence 
viscosity can also be destabilizing. 

An argument due to C. C. Lin (MIT) is sketched below. 

Near the wall of a parallel shear flow such as Blasius boundary layer of thickness S (or a 
channel flow of depth 5), a wave- like disturbance introduces an oscillatory Stokes boundary 
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Figure 5.4.1: An oscillatory Stokes boundary layer 

layer. Refering to Figure 5.4.1, the oscillatory boundary layer has the thickness of the order 
6o = ^J2u/iJJ. Let us calculate the Reynolds stress due to wave disturbances in this boundary 
layer. 

For this purpose let us consider first the vorticity disturbance in the Stokes boundary 
layer. After linearization, the vorticity satisfies 

C; + UC + v'Cy = + Cyy) (5-4.5) 

In the case of traveling wave disturbance 

Inside the Stokes layer, C,' = —u'y so that u' ~ 5oC,' ■ From continuity we have v' ~ fc^oit'. Now 
consider the two convection terms: 

ukc u 

CI ojC c 

where C — oj/k, and 

Close enough to the wall U/C 1; the convective terms are negligible. 
On the other hand, 

Cxx/Cyy ~ {kSf 

Hence in the wall boundary layer, we have approximately 

C = <y (5.4.6) 

The boundary conditions are : 

l^' ^ Ae'^kx-u,t) ^ y = (5.4.7) 
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and 



(5.4.8) 



This boundary value problem is just the Stokes problem of a plate oscillating in its own 
plane. The solution is, in complex form, 



C' = [Ae'^e'^y) , where 9 ^ kx - cot, a = 



(5.4.9) 



The real part is to be taken at the end. 
Let us calculate u' from the vorticity. 
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and v' from continuity, 

v'^- [\'dy = -—e''^ [\l - e^y)dy = - 
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For y/^o <C 1 or ay <C 1 we can approximate n', v' as power series 

A -n 
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(5.4.10) 



(5.4.11) 



(5.4.12) 



(5.4.13) 



Now we shall take the time average to get the Reynolds Stress u'v' . For two sinusoidal 
functions / = 3?Fe~*'^* and g — SRGe*"*, it can be shown that the time average over a period 
is 
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t+27r/w 1 1 

fgdt = -^{FG*) = -^{F*G) 



(5.4.14) 



Here the time average is equal to the space average over a wavelength and to the phase 
average. 

Substitituting (5.4.12) and (5.4.13) into (5.4.14), we get 
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-u'v' = --3? 
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The Reynold stress is therefore positive 

-pu'v' = ^-TTT^yS > (5.4.15) 
24 

to the leading order and is inversely proportional to On the other hand, the rate of 
dissipation is proportional to the integral of 

fiC'^ oc ii\A\^ (5.4.16) 

For sufficiently small viscotiy work done by the Reynolds stress to the mean shear can 
therefore overwhelm disspation, leading to instatbilty. This estimate has also been cofirmed 
by fuller anslysis (see Lin, p 62ff.) . 

Homework Use the solution (5.4.9), (5.4.10) and (5.4.11) without assuming small y/(5o 
and compute pM and iiN and discuss the result. 



5.4.2 Orr-Sommerfeld theory 

The comoplete solution of the viscous instability problem begins with the linearized Navier- 
Stokes equations. Starting from a basic shear flow Q = {U (y), 0), we consider an infinitesimal 
disturbance with the velocity field q — {u{x,y,t),v{x,y,t)) which is much smaller than U. 
Introduce the stream fundtion 

dip dtp 
dy ' dx 

and assume a normal mode disturbance 

^ = /(y)e*'=(^-^*). where C = u/k 
subject to the boundary conditions that 

^^CV'y^O, y = 0, onthewah (5.4.17) 

^ = 0, ^/-y = 0, y = d{or oo) (5.4.18) 

After nondimensionalization and lineaization of the Navier-Stokes equations we get the 
Orr-Sommerfeld equation for f{y) 

f" _ 2/cV'" + k^f = -ikR [{U{y) - C){f" - k^f) - U"f] (5.4.19) 
where R is the Reynolds number. The boundary conditions are 

/ = /' = 0, y = 0,l (5.4.20) 

the eignevalue problem can be solved numerically. 

For the boundary layer, computations (Tolmien and Schlicliting) show that the threshold 
of instability C, = is a hairpin-like curve in the plane of kSi{x) vs. Ri{x) = U6i/u where 
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6 - displacement thickness 
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Figure 5.4.2: Instbility region of a Blasius boundary layer, fc^^i vs. Re^^ where S\ is the 
displacemenet thickness which increases with J {x) . 




5i(,x) oc ^/x is the displacement thickness, as shown in Figure 5.4.2. At a fixed x, a sinudoidal 
disturbance of wave number k is stable if the Reynolds number is within certain hairpin-like 
band. The is also a critical Reynolds number below which all disturbaqnces are stable, and 
above which some distrubances are unstable. 

Verifications were first made in a landmark experiment by Schubauer, Skramstad and 
Dryden (1945). With controled disturbances where the amplitude and wavelength of the 
sinusodial distrubance were introduced by a loudspeaker, the Tolmien-Schlichting threshold 
was confirmed. Under the test condition 



they found that natural disturbances were still wave-like at a; = but waves are amplified 
at 9 ft. Occasional bursts appeared at 10.5 feet and signals become random at 11 ft. Later 
experiments and nonlinear theories showed the following development as 
X increases : 

1. For small x a laminar boundary layer is stable. 

2. For sufficiently long waves, linear instability leads to Tolmien-Schlichting waves. 

3. Tolmien-Schlichting waves grow large. Nonlinearity sets in. Disturbances become 3- 
dimensional. Spanwise periodicity appears. Local instablitiy occurs. 
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4. Reynolds-stress generates higher harmonics. Turbulent spots appear. 




Figure 5.4.3: Comparison of theory and experiments on boundary layer instability. 

5. Turbulent spots grow and overlap. 

6. Boundary layer becomes turbulent. 

This evolution is summarized in Figure 5.4.5 
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Figure 5.4.5: Spatial evolution of Blasius' boundary layer from instability to turbulence 
instability. 



